Dynamical structure factors of the magnetization-plateau state in the 5 = 1/2 
bond-alternating spin chain with a next-nearest-neighbor interaction 
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We calculate the dynamical structure factors of the magnetization-plateau state in the S = 1/2 
bond-alternating spin chain with a next-nearest-neighbor interaction. The results show characteristic 
behaviors depending on the next-nearest-neighbor interaction a and the bond-alternation 5. We 
discuss the lower excited states in comparison with the exact excitation spectrums of an effective 
Hamiltonian. From the finite size effects, characteristics of the lowest excited states are investigated. 
The dispersionless mode of the lowest excitation appears in adequate sets of a and 5, indicating 
that the lowest excitation is localized spatially and forms an isolated mode below the excitation 
continuum. We further calculate the static structure factors. The largest intensity is located at 
q = n for small S in fixed a. With increasing S, the wavenumber of the largest intensity shifts 
towards q — it/2, taking the incommensurate value. 

PACS numbers: 75.25.-l-z, 78.70.Nz, 75.10.Jm 



I. INTRODUCTION 

One-dimensional (ID) quantum spin systems have at- 
tracted great attention both theoretically and experimen- 
tally. Fascinating aspects can be seen in some char- 
acteristic systems under magnetic fields. Recently, it 
was shown numerically that there appears a plateau re- 
gion on the magnetization curve at half of the satura- 
tion value in the S = 1/2 bond-alternating spin chain 
with a next-nearest-neighbor (NNN) interaction 
In this half-magnetization-plateau region, furthermore, 
the wave function of the lowest energy state was shown 
to have the periodicity four in units of the lattice con- 
stant, although the periodicity of the original Hamilto- 
nian is two 0, H, 1). The result indicates that sponta- 
neous translational-symmetry breaking occurs ^. 
This feature is consistent with the necessary condition 
for the appearance of the half-magnetization-plateau in 
the S = 1/2 spin chain [||. The stabihty of the 
half-magnetization-plateau was discussed using bosoniza- 
tion technique combined with the renormalization-group 
theory and the perturbation calculation On the 

basis of the results, the phase diagram of the half- 
magnetization-plateau, which is consistent with the nu- 
merical results 1^, 1^ , was obtained as a function of the 
parameters of the bond-alternation and the NNN inter- 
action [||. A simple picture of the half-magnetization- 
plateau state was presented, indicating the two-fold de- 
generate states with singlet and triplet pairs occupying 
strong bonds alternately Furthermore, it was shown 
that the elementary excitation is given by a superposition 
of massive kinks and anti-kinks which have — ±1/2, 
respectively 

The dynamical structure factor (DSF) provides us 
with important information about the intensity of the 
magnetic excitation as a function of energy-momentum 
transfer. Since the DSF can be observed experimen- 
tally by inelastic neutron-scattering measurements, a de- 
tailed comparison between theoretical and experimental 



results has been developed. The DSF of the S = 1/2 
bond-alternating spin chain with a NNN interaction in 
magnetic fields was investigated ^, ^ 0| in connection 
with experimental findings of the spin-Peierls material 
CuGeOa ^ , since this model with adequate parameters 
can be an effective model for CuGeOs 11^. Critical 
and dynamical properties of typical (quasi-jlD S = 1/2 
spin-gapped systems such as the bond-alternating spin 
chain, the spin chain with a NNN interaction, and the 
two-leg ladder system were investigated systematically 
on the basis of bosonization technique 1 11 . It was shown 



that the DSF exhibits different peak structures among 
these three systems, when the gap is collapsed by mag- 
netic fields. The DSF of the S = 1/2 spin chain with 
a NNN interaction was calculated at half of the satu- 
rated magnetization using the numerical diagonalization 
method The distribution of the intensity was dis- 

cussed, putting stress on the behavior around the field- 
induced soft mode. One of the present authors calculated 
the DSF of the half-magnetization-plateau state in the 
5 = 1/2 bond-alternating spin chain with a NNN inter- 
action using two sets of parameters It was shown 
that the periodicity and the amplitude of the dispersion 
relation of the transverse DSF are sensitive to the pa- 
rameters. In spite of this investigation, the dynamical 
property of the half-magnetization-plateau state has not 
been fully understood. A detailed investigation is desir- 
able. 

In this paper, we investigate dynamical properties of 
the half-magnetization-plateau state in the S = 1/2 
bond-alternating spin chain with a NNN interaction. Ac- 
cording to the phase diagram ||, ^, ||, we calculate the 
DSF systematically in the parameter space of the half- 
magnetization-plateau. In §2, we briefly summarize the 
method for the numerical calculation. In §3, we show the 
results for the transverse DSF. In the half-magnetization- 
plateau state with a strong bond-alternation and a weak 
NNN interaction, the Hamiltonian of the system can be 
mapped onto an effective Hamiltonian described by the 
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ID S = 1/2 Heisenberg-Ising model in zero magnetic 
field m . We compare the numerical results with the 
exact excitation spectrums of the effective Hamiltonian. 
We also estimate the finite-size effects carefuhy to in- 
vestigate characteristics of the lowest excited states. We 
further investigate the static structure factor (SSF) in the 
half-magnetization-plateau state. Section IV is devoted 
to the summary of the paper. 



II. MODEL AND METHOD 

We consider the S" = 1/2 bond-alternating spin chain 
with a NNN interaction in magnetic fields described by 
the following Hamiltonian, 



H 

Ho 



Ho+TLz, (1) 

N 



N 



Hz = -gtXBHY^St 



(3) 



1=1 



where 5 denotes the bond-alternation, a denotes the 
NNN interaction, N is the total number of the site, and 
H is magnetic field. We set J = 1 and gfj^B — 1- The 
periodic boundary condition is applied. 

In magnetic field along the z axis, rotational symme- 
try around the x and y axes is broken, while that around 
the z axis remains. Therefore, the Hamiltonian can be 
classified into the subspace according to the magnetiza- 
tion m = M/N with M = J2f=i We set that l^-o,,™) 
is the eigcnfunction of the lowest eigenvalue of Ti. with 
magnetization m. In given m, the DSF can be expressed 
in the form of the continued fraction [[13 as 
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{^O,m\SfS^\^0,,n), 
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where 5*^ 



w-l-Je-l-i^^ii/, h. — \, and is the lowest eigenvalue of 
Tig with magnetization m. In Eq. {^, ai and (i — 
0, 1, 2, • • •) are the diagonal and the sub-diagonal elements 
of the tridiagonalized Hamiltonian TYq. We put the lat- 
tice constant to unity. The DSF S'^^q^iS) is obtained by 
the relation S'^(g,w) = [5'+(g,w) -f S'~(q, w)] /4, where 



S'^{q,u!) = Sy{q,u!) ^ S^{q,uj). We can determine poles 
and residues of this continued fraction numerically. The 
total contribution of C'^{q, uj) for fixed q is normalized to 
unity, because the following sum rule has to be satisfied. 



S^iq) 



dujS''{q,uj). 



(7) 



To calculate S''^(g,a;) according to (Q), we have to set 
the magnetic field Hp of the half-magnetization-plateau 
region. We now express as (TV, M) in the A^-spin 
system, fixing m = 1/4. We define the magnetic field 
where the magnetization changes from M to M± 1 as ll^ 



H+ = E°{N, M -f 1) - E°{N, M), 
= E°{N, M) - £;"(iV, M - 1), 



(8) 
(9) 



where and H~ are the magnetic fields at the high- 
est and the lowest ends of the magnetization plateau, 
respectively [|], |l^. Using 7J+ and we can express 
the magnetic field Hp. Here, we set Hp = +H~]/2. 

We first calculate |*o,m), E^{N,M), and E^{N,M±1) 
numerically with M — N /A, fixing TV. The Hamiltonian 
is then tridiagonalized by the Lanczos method. We sub- 
stitute Hp, ai and for the expression (6). Instead of 
taking the limitation e — > -1-0, we set e = 3.0 x 10~^. The 
convergency of the continued fraction can be estimated 
quantitatively by the modified Lenz method ||l8|. The 
finite-size effects are estimated in the same way as that 
used in Refs. 18 and 19. In this way, we obtain the DSF 
in the half-magnetization-plateau state. 



III. NUMERICAL RESULTS 

A. DSF in the half-magnetization-plateau state 

In Fig. 1, the transverse DSF S^{q,uo) in the half- 
magnetization-plateau state is shown for N=28 and M — 
7. The intensity of S^{q, lo) is represented using the same 
scale, and is proportional to the area of the full circle. 
The accuracy of the continued fraction in (6) becomes 
worse, when we treat the higher energy region. To get 
accurate results, we have increased the order of the con- 
^binned fraction. As a result, the continued fraction has 

— converged within the relative error 0(10"^") for lo < 1, 

-and 0(10-2) for w > 3. 

_ In larger 5 and/or smaller a, the separation of the ex- 
■ ' citation bands becomes apparent, and we can see nearly 
(Qihree excitation bands in the region w < 1, w ~ 2, and 
ti> > 3. The width of the excitation band strongly de- 
pends on the parameters a and 6. Furthermore, narrow 
and almost dispersionless bands appear in some adequate 
sets of a and 5 close to 2q: + 5 = 1 . Concerning this fea- 
ture, we will show details later. The larger intensity lies 
mainly in the lowest excitation bands. The wavenumber 
of the largest intensity in given a shifts from tt to 7r/2, 
as 5 increases. The results indicate that the effect of the 
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bond-alternation becomes dominant. Thus, we turn our 
attention to the behavior of the lowest excitation band. 

We have also calculated the longitudinal DSF S^{q,uj), 
and have found that S^{q,uj) does not exhibit noticeable 
features depending on the parameters a and S. There- 
fore, we consider only the behavior of in the 
following. 



B. Comparison with the excitation spectrums of 
the effective Hamiltonian 

We first compare the distribution of the intensity of 
S^{q,u!) in the lowest excitation band with the exact ex- 
citation spectrums for the effective Hamiltonian. In the 
half-magnetization-plateau state with 5 ~ 1 and a ^ 0, 
the model described by (0) can be mapped onto the ID 
S = 1/2 Heisenberg-Ising model in zero field described 
as i,|l7 



^Jxy{SjSj^i + S^SJ_^_^) + JzSjSj_^_i^ , (10) 



where Xy = [2a - (1 - 5)]/2 and ^ [2a + (1 - 5)]/ A. 
The effective spin- 1/2 operators are S'f = 2[S2i_i — 
1) = 2(5|, - 1),^+ = V2S+_, - ~^/2S+, and Sr = 

y/2S~i_^ = -V2S'2l, {I = 1,2,3, •••,iV/2). The ground 
state of ([l0| ) is a Neel-like state with two-fold degener- 
acy, which corresponds to the singlet and triplet pairs 
occupying strong (1-1- (5)-bonds alternately in the orig- 
inal Hamiltonian [^|. The elementary excitation, which 
keeps = ±.1 and 0, can be given by a superposition 
of massive kinks and anti-kinks which have S*^ = ±1/2, 
respectively. 

In Fig. 2, the results for N=28 and M = 7 are shown. 
The solid lines represent the exact bounds of the ele- 
mentary excitation for the effective Hamiltonian obtained 
by the Bethe ansatz method |21, In smaller 5 of 

given a, the intensity distributes beyond the upper and 
lower bounds of the excitation continuum for the effec- 
tive Hamiltonian. In larger 5 of given a, on the other 
hand, the intensity distributes within the excitation con- 
tinuum for the effective Hamiltonian. In the latter case, 
the ground state can be given by the singlet and triplet 
pairs occupying strong (1 -I- (5)-bonds alternately, and the 
elementary excitation can be described by a superpo- 
sition of massive two kinks and anti-kinks which have 
= ±1/2. The schematic picture for the elementary 
excitation is the same as that shown in Fig. 5 of Ref. 3, 
where the kink and anti-kink are represented by domain- 
wall-like excitations from the ordered ground state men- 
tioned above. In the former case, it is not easy to present 
a schematic picture for the elementary excitation. The 
higher order terms which are released to obtain the effec- 
tive Hamiltonian may become important. It is interesting 
that even in a > 0.1, the intensity distributes within the 
excitation continuum for the effective Hamiltonian, al- 



though the mapping onto the effective Hamiltonian can 
be successful for 5 ~ 1 and a ~ 0. 

In the parameter region where 2a + 6 = 1 is satisfied, 
the effective Hamiltonian is reduced to the ID 5 = 1/2 
Ising model, and then the elementary excitation shows 
a dispersionless mode. Note that this parameter region 
in zero field is called the Shastry-Sutherland line P3[ . 
In Fig. 2, we can see S^{q,uj) in the parameter region 
2a + 6 ^ 1; {a,S) = (0.1, 0.8), (0.2, 0.6), (0.3, 0.4) and 
(0.4,0.2), where the dispersionless mode of the effective 
Hamiltonian is represented by the broken line in given a. 
At S = 0.8 and a = 0.1, the agreement between the nu- 
merical results and the results from the effective Hamil- 
tonian is excellent, indicating that the lowest excitation 
is described by a superposition of massive two kinks and 
anti-kinks, and is well localized in the real space. As the 
parameter region is removed from the successful mapping 
point where (5 ~ 1 and a ~ 0, the intensity spreads more 
widely and the spatial localization of the elementary ex- 
citation becomes worse. 

To see the characteristics of the low-lying excitations 
in detail, we next investigate the finite-size effects of the 
lowest and the next-lowest excited states. According to 
Refs. 18 and 19, we investigate the size dependence of 
the poles and their residues of the continued fraction 
C^{q,uj), to discuss whether a pole belongs to an ex- 
citation continuum or forms an isolated mode. When a 
pole belongs to an excitation continuum, at least either 
of its position or its residue has appreciable size depen- 
dence. On the other hand, when a pole forms an isolated 
mode, the position and the residue of a pole has little 
size dependence. This method was successfully used to 
investigate the isolated mode of the S = 1 Haldane gap 
system Later on, characteristics of the lowest ex- 

cited state for the S" = 1/2 bond-alternating spin chain 
with a NNN interaction was investigated using the same 
method |2^. In this study, it was shown that a defi- 
nite conclusion concerning the isolated mode is difficult 
to draw only from this procedure. 

In Fig. 3, we show the typical results for the system- 
size dependence of poles and residues of the contin- 
ued fraction in the present system. Note that we have 
confirmed that the same system-size effects already ob- 
tained for the poles and residues of the 5=1/2 XXX 
chain jl^, |2^ can be reproduced. 

We first consider the behaviors in the region where the 
intensity spreads within the excitation continuum of the 
effective Hamiltonian. The typical results are shown in 
Fig. 3 (a), which corresponds to the bottom figure in Fig. 
2 (c). In this region, the size dependence of the residues 
is not noticeable. The positions of the poles for the lowest 
excited states show almost no size dependence. Yet, the 
positions of the poles for the next lowest states decrease, 
as the system size increases. Judging from these results 
and the fact that the intensity spreads within the excita- 
tion continuum, the lowest mode of the present region is 
the lower bound of the excitation continuum. 

We next consider the behaviors in the region where the 
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intensity spreads beyond the excitation continuum. In 
this region, it is difficult to present a definite conclusion 
whether the lowest excited states form an isolated mode. 
As shown in Figs. 3 (b) and (c), which correspond to the 
top figures in Figs. 2 (b) and (d) respectively, we cannot 
see clear size dependence of the residues. The positions 
of the poles for the lowest excited states scarcely show 
size dependence. Nevertheless, the positions of the poles 
for the next lowest states decrease around q = 0, 7r/2, 
and TT, or around two wavenumbers among them, as the 
system size increases. Thus, in this parameter region, 
the lowest excited states may become a lower bound of 
the excitation continuum around q = 0, 7r/2, and tt, or 
around two wavenumbers among them, while they may 
form an isolated mode in the other wavenumbers. 

In the above arguments, we leave the results for the 
parameter region where the almost dispersionless mode 
appears. Analysis in such parameter regions will be per- 
formed in the next subsection. 



C. DSF in the dispersionless mode 

In the parameter region away from 5 1 and a ~ 
0, the almost dispersionless mode of S''^{q,uj) appears 
around 6 ~ 0.64 in a = 0.2, 6 ^ 0.5 in a = 0.3, and 
6 ~ 0.38 in a = 0.4 as shown in Fig. 4. The size de- 
pendence of poles and residues at these points is shown 
in Fig. 5. Note that the similar behaviors are seen be- 
tween the cases of {a,S) — (0.1,0.8) and (0.2,0.64), and 
between {a, 6) — (0.3,0.5) and (0.4,0.38), respectively. 
The residues in {a, 6) = (0.1,0.8) and (0.2,0.64) show 
no size dependence, while those in (a, 6) = (0.3, 0.5) and 
(0.4, 0.38) decrease around q ^ 0.757r and n with increas- 
ing the system size. The positions of the poles in (a, S) = 
(0.1, 0.8) and (0.2, 0.64) show no size dependence. There- 
fore, in {a, 5) = (0.1,0.8) and (0.2,0.64), the lowest ex- 
citation is localized spatially and forms apparently the 
isolated mode below the excitation continuum. In fact, 
no intensity distributes between the lowest dispersionless 
mode and the excitation band around cu ~ 2.0 in the fig- 
ure of {a, 5) = (0.1,0.8) in Fig. 1 (a). Note that the 
same occurrence is seen also in (a, <5) = (0.2,0.64), al- 
though we have not presented the results in Fig. 1 (b). 
On the contrary, in {a,S) = (0.3,0.5) and (0.4,0.38), the 
positions of the poles for the lowest excited states show 
no size dependence, while those for the next lowest states 
decrease around q ~ 0.757r and vr with increasing the sys- 
tem size. Judging from the results, the lowest excitation 
in (a, S) = (0.3, 0.5) and (0.4, 0.38) is probably the lower 
edge of the excitation continuum around q ~ 0.757r and 
TT, and is the isolated mode in the other wavenumbers. 

In this study, we have neither found the completely 
isolated mode in a = 0.3 nor 0.4. Nevertheless, there 
may appear the completely isolated mode close to the 
parameters (a, 5) = (0.3,0.5) and (0.4,0.38). We leave 
this issue in the future. 

In the previous and present subsections, we have dis- 



cussed the isolated mode between the ground state and 
the excitation continuum. The isolated mode can be seen 
in several ID spin-gapped systems such as the Haldane- 
gap system [^9[ gj, |2^, the S = 1/2 bond-alternating 
spin chain with a NNN interaction [b| BO i |2^ , P7| , PSj , the 
5 = 1/2 two- leg spin ladder systenr|29|, pO|, |31|, |32||, and 
so on In some systems, the isolated mode was 

observed experimentally |2^, 0. It may be inter- 
esting that such an isolated mode has appeared also in 
the field-induced ordered state of the ID quantum spin 
system. 



D. SSF in the half- magnetization-plateau state 

We calculate the SSF S^{q) in the half-magnetization- 
plateau state using the sum rule (|^). The results are 
shown in Fig. 6. The finite-size effects appear only at 
q = n OT q = tt/2 where the intensity takes the maximum. 
The largest intensity is located at g = tt in the small d 
region. As d increases, the wavenumber of the largest 
intensity shifts towards q = tt/2, taking the incommen- 
surate value of q. This occurrence reflects that the bond- 
alternation becomes effective. In a = 0.4, the region of 
the incommensurate wavenumber becomes wider as com- 
pared to that for other a's. 

In other ID spin-gapped systems such as the Haldane- 
gap system |l^ and the S* = 1/2 Heisenberg chain with 
a NNN interaction |Q, the incommensurate mode be- 
tween 7r/2 and tt is difficult to observe. The dominant 
incommensurate mode of the transverse spin correlation 
may be characteristic of the magnetization-plateau state 
in the present system. 



IV. SUMMARY 

We have investigate the DSF of the half-magnctization- 
plateau state in the S* = 1/2 bond- alternating spin chain 
with a NNN interaction. The results are summarized in 
the phase diagram of Fig. 7. The half-magnetization- 
plateau appears along the thick lines in given a [Q . At 
the full circles, the intensity distributes within the excita- 
tion continuum for the effective Hamiltonian, while at the 
open circles, the intensity distributes beyond the excita- 
tion continuum. In the former regions, the ground state 
and the elementary excitation are probably described 
by a simple picture based on the effective Hamiltonian. 
From the analysis of the size dependence, we have con- 
cluded that the lowest excited states in the former region 
are the lower bound of the excitation continuum, while 
those in the latter region may form the isolated mode be- 
low the excitation continuum. Along the dotted line, the 
dispersionless behavior comes out, indicating that the el- 
ementary excitation is localized spatially. From the size 
dependence, we have concluded that the dispersionless 
mode is isolated below the excitation continuum. 
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Above the upper solid line, the largest intensity of SSF 
is located at q = tt/2, and below the lower solid line, the 
largest intensity of SSF is located at q = n. As 6 in- 
creases between two solid lines, the wavcnunibcr of the 
largest intensity shifts from q — tt to q — 7r/2, taking the 
incommensurate value. The results indicate that with in- 
creasing S in fixed a, the dominant mode of the transverse 
spin correlation shifts tt to 7r/2, taking an incommensu- 
rate value. 



sions. Our computational programs are based on TIT- 
PACK Ver. 2 by H. Nishimori. Numerical computation 
was partly carried out at the Yukawa Institute Computer 
Facility, Kyoto University, and the Supercomputer Cen- 
ter, the Institute for Solid State Physics, University of 
Tokyo. 



V. ACKNOWLEDGMENTS 



We would like to thank I. Harada, M. Takahashi and 
T. Tonegawa for useful comments and valuable discus- 



T. Tonegawa, T. Nishida, and M. Kaburagi, Physica B 
246-247 (1998) 368; private communication. [19 
T. Tonegawa, T. Hikihara, K. Okamoto, and M. [20 
Kaburagi, Physica B 294-295, 39 (2001). 

K. Totsuka, Phys. Rev. B 57, 3454 (1998). [21 
M. Oshikawa, M. Yamanaka, and I. Affleck, Phys. Rev. 
Lett. 78, 1984 (1997). [22 

D. Poilblanc, J, Riera, C. A. Hayward C. Berthier, and 
M. Horvatic, Phys. Rev. B 55, R11941 (1997). [23 
S, Mori and I. Harada, J. Phys. Chem. Solids 60, 1133 

(1999) . [24 
W. Yu and S. Mass, Phys. Rev. B 62, 344 (2000). 
Y. Fagot-Revurat, M. Horvatic, C. Berthier, J. [25 
P. Boucher, P. Segransan, G. Dhalenne, and A. 
Revcolevschi, Phys. Rev. B 55, 2964 (1997). [26 
J. Riera and A. Dobry, Phys. Rev. B 51 16098 (1995). 
G. Castilla, S. Chakravarty, and V. J. Emery, Phys. Rev. [27 
Lett. 75 1823 (1995). 

R. Chitra and T. Giamarchi, Phys. Rev. B 55, 5816 
(1997). ' [28; 

C. Gerhardt, A, Fledderjohann, E. Aysal. K.-H. Miitter, 
J. F. Audet, and H. Krogner, J. Phys.: Condens. Matter [29 
9, 3435 (1997). 

M. Usami and S. Suga, J. Phys. Soc. Jan. 69, 3419 

(2000) . ' [30 
F. Mila, Eur. Phys. J. B 6, 201 (1998). 

E. R. Gagliano and C. A. Balseiro, Phys. Rev. Lett. 59, [31 
2999 (1987). 

[16] To calculate the precise values of the magnetic fields in 
the thermodynamic limit at both ends of the magnetiza- 
tion plateau, finite-size-scaling analysis based on confor- [32 
mal field theory is worthwhile as shown in Refs. |l| p^ . 
Nevertheless, we have neglected the finite-size correction [33 
in ^ and (j^, because the error in the magnetic fields [34 
thus obtained shifts the DSF slightly along the ui axis. 

[17] T. Sakai and M. Takahashi, Phys. Rev. B 43, 13383 
(1991); J. Phys. Soc. Jpn. 60, 3615 (1991). 

[18] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. 
P. Flannery, Numerical Recipes, (Cambridge University 



Press, 1992). 

M. Takahashi, Phys. Rev. B 50, 3045 (1994). 

H. Yokoyama and Y. Saiga, J. Phys. Soc. Jpn. 66, 3617 
(1997). 

A. H. Bougourzi, M. Karbach, and G. Miiller, Phys. Rev. 
B 57, 11429 (1998). 

M. Takahashi, Thermodynamics of One-Dimensional 
Solvable Models, (Cambridge University Press, 1999). 

B. S. Shastry and B. Sutherland, Phys. Rev. Lett. 47, 
964 (1981). 

S. Ma, C. Broholm, D. H. Reich, B. J. Sternlieb, and R. 
W. Erwin, Phys. Rev. Lett. 69, 3571 (1992). 

I. A. Zaiiznyak, S.-H. Lee and S. V. Petrov, Phys. Rev. 
Lett. 87, 017202 (2001). 

G. S. Uhrig and H. J. Schuiz, Phys. Rev. B 54, R9264 
(1996). 

M. Ain, J. E. Lorenzo, L. P. Regnault, G. Dhalenne, A. 
Revcolevschi, B. Hennion, and Th. Jolicoeur, Phys. Rev. 
Lett. 78, 1560 (1997). 

P. V. Shevchenko, V. N. Kotov, and O. P. Sushkov, Phys. 
Rev. B 60, 3305 (1999). 

O. P. Sushkov and V. N. Kotov, Phys. Rev. Lett. 81, 
1941 (1998); V. N. Kotov, O. P. Sushkov and R. Eder, 
Phys. Rev. B 59, 6266 (1999). 

C. Jurecka and W. Brenig, Phys. Rev. B 61, 14307 
(2000). 

M. Windt, M. Grueninger, T. Nunner, C. Knetter, K. 
Schmidt, G. S. Uhrig, T. Kopp, A. Freimuth, U. Am- 
merahl, B. Buechner, and A. Revcolevschi, Phys. Rev. 
Lett. 87 127002 (2001). 

C. Knetter, K. P. S chmidt, M. Griininger, and G. Uhrig, 
cond-mat/0106077[ 

K. Damle and S. Sachdev, Phys. Rev. B 57, 8307 (1998). 
S. Trebst, H. Monien, C. J. Hamer, Z. Weihong, and R. 
R. P. Singh, Phys. Rev. Lett. 85, 4373 (2000); W. Zheng, 
C. J. Hamer. R. R . P. Singh, S. Trebst, and H. Monien, 
cond-mat/0010354 



This figure "figlab.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0110284vl 



This figure "figlcd.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0110284vl 



This figure "fig2ab.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0110284vl 



This figure "fig2cd.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0110284vl 



This figure "fig3.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0110284vl 



This figure "fig4.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0110284vl 



This figure "fig5.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0110284vl 



This figure "fig6.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0110284vl 



This figure "fig7.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0110284vl 



